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Cu NQR Evidence for Spatial Variation of Hole Concentration in La2_2^^Sra,Cu04 
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We report experimental evidence for the spatial variation of hole concentration Xhoie in the high 
Tc superconductor La2-i:Sra;Cu04 (0.04 < a:: < 0.f6) by using ''^Cu NQR for ''^Cu isotope enriched 
samples. We demonstrate that the extent of the spatial variation of the local hole concentration 
Axhoie is reflected on ^^l/Ti and deduce the temperature dependence. Axhoie increases below 
500 — 600K, and reaches values as large as Axhoie/x 0.5 below 150K. We estimate the length 
scale of the spatial variation in x^oie to be Rhoie ^ 3 nm from analysis of the NQR spectrum. 
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The mechanism of high Tc superconductivity remains 
enigmatic even after 15 years of its discovery [Q. It is well 
known that the central idea in the search of high ma- 
terials was to look for materials with polaronic effects 
caused by coupling between electrons and the lattice. 
However, most recent theoretical debates are based on 
the assumption that the holes are homogeneously doped 
into Cu02 planes, even if hole-doping is achieved by the 
substitution of ions with different ionicity which in effect 
creates alloys. A clear and widely accepted counter exam- 
ple against such a simplified picture is the stripe phase in 
Nd co-doped La2-2;Sr2;Cu04, where the Coulomb poten- 
tial from the distorted lattice slows down spin and charge 
density waves |^ . It is also known that the high mobility 
of the oxygen atoms results in electronic phase separation 
[^-^ between the superconducting and antiferromagnetic 
phase in super-oxygenated La2Cu04+5 as evidenced by 
^^^La NMR (nuclear magnetic resonance) [||. A natural 
question to ask is whether similar situations arise in other 
materials albeit in a less robust manner 0. In fact, there 
is ample evidence for some sort of spatial inhomogeniety 
in the Cu02 plane from earlier NMR and NQR (nuclear 
quadrupole resonance) studies in La2-xSr^Cu04 
such as splitting of the ^^Cu NQR and NMR lines due 
to ineqivalent Cu sites with different EFG (electric field 
gradient) tensors H,^, drastic broadening of the ^^Cu 
NQR line at low temperatures [||jll), and ^^Cu NMR Hne 
broadening due to the modulation of orbital shifts JTif . 
More recent studies using ^'^Cu NQR and NMR wipe- 
out have characterised the glassy nature of the 
slowing down of the stripe inhomogeneity below tempera- 
tures ~100K. However, no clear picture has emerged that 
discerns and relates the effects of genuine electronic phase 
separation, stripe modulation, and random substitution 
of donor ions. On the other hand, a recent STM (scan- 
ning tunnelling microscopy) study on the surface state of 
Bi2Sr2CaCu208+5 cleaved at low temperatures revealed 
nano-scale spatial variations of the electronic state [^ . 
Whether such nano-scale modulations are universally ob- 
servable in the bulk and other high Tc materials remains 
to be seen, but the STM results have enhanced the inter- 



est in the potential impact of the spatial inhomogeniety 
of the electronic properties in the Cu02 plane. 

In this Letter, we utilize a new simple trick based on 
the local nature of the ^'^Cu NQR technique to probe the 
spatial inhomogeniety of the Cu02 plane. Unlike scatter- 
ing techniques, NQR does not require spatial coherence 
over tens of nm's, hence it is an ideal technique to probe 
spatial variations in the local hole concentration Xhoie at 
short length scales. Our main results are presented in 
Fig. 1, where we plot the temperature dependence of the 
extent of the spatial variation Axhoie as measured by 
the distribution in (data with solid line) and by 

analysis of the NQR spectrum (dashed lines). Contrary 
to common assumptions that doped holes are uniformily 
distributed in the Cu02 planes (i.e. Aa;;io/e=0), we find 
that a non-zero value of Axhoie exists at all tempera- 
tures and furthermore shows a factor 2-3 increase from 
500-600K down to ~150K. We also demonstrate that the 
NQR spectrum is consistent with a model where the spa- 
tial variation in Xhoie takes the form of patches in the 
Cu02 plane, where some patches are more metallic and 
some more insulating, with a patch radius Rhoie nm. 

We shall now describe the details of the proceedure 
used to obtain Ax hole shown in Fig. 1. The nuclear spin- 
lattice relaxation rate is given by the formula: 
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where w„ is the NQR frequency, ^^A{q) is the hyperfine 
form factor, and ^(q, a;„) is the dynamic stucture fac- 
tor of the Cu electron spins measures the 
local low frequency Cu spin fluctuations at w^. Because 
of the experimental ease, the temperature dependence of 
is generally measured at the peak of the spectrum. 
However, in various materials including La2-a:Sra:Cu04 
pot , 1/Ti often depends on frequency across the reso- 
nance spectrum due to a spatial distribution of the elec- 
tronic states. A typical example is shown in Fig. 2a for 
X = 0.115. For convenience we have defined l/Tf, 1/T^, 
l/T{' and 1/T^^^^^ as ^^1/Ti measured at the peak, up- 
per half intensity, lower half intensity and lower one tenth 
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intensity of the A-line respectively. 

The temperature dependence of l/Tj^, and 
llj^-^/w piottejj Yig. 3. All data were taken in 
the temperature region above Cu wipeout , where full 
NQR signal intensity is observable. For comparison, Fig. 
3 also shows 1 /T" for a variety of samples which all show 
agreement with previous work We note that all of 
our results are insensitive to RF pulse width. The most 
surprising discovery of the present work is that 
l/T]^ and ^^^"^ show qualitatively different temper- 

ature dependence. For example, in the x = 0.07 sample, 
^^l/Ti measured at the half intensity of the lower (up- 
per) frequency side of the NQR spectrum exhibits semi- 
quantitatively the same behaviour as for x = 0.04 

{x — 0.115). This is consistent with the fact that the 
lower (upper) frequency side of the NQR spectrum for 
X = 0.07 roughly coincides with the peak NQR frequency 
of X = 0.04 (a; = 0.115) as shown in Fig. 2b. Our finding 
immediately establishes that within a single sample with 
a fixed nominal hole concentration x, there are some seg- 
ments with higher and lower hole concentrations Xhoie- 

By comparing 1 /T^ with smoothly interpolated values 
of 1/Ti for various samples, we can estimate the devia- 
tion t\xhoie of local hole concentration from the spatial 
average, x. For example, l/Tj^ for x = 0.07 at 200K is 
close to l/Ti* for x ~ 0.045, implying that the lower fre- 
quency side of the x — 0.07 spectrum corresponds to res- 
onance from segments with l\xhoie —0.025. Similarly, 
l/r+ for X = 0.07 is close to l/Tf for x ~ 0.10, implying 
the upper frequency side corresponds to segments with 
^Xfioie ~ -1-0.03. Using instead results in an 

overall 30 — 40% increase in Axhoie and does not effect 
our conclusions. Another important feature in Fig. 3 is 
that cuts through the lines of with decreas- 

ing temperature. This implies that Axhoie increases with 
decreasing temperature, as shown in Fig. 1. Since the 
magnitude of Axhoie estimated from and 1/Tf is 

identical within experimental uncertainties, we plot the 
magnitude of Axhoie ■ To the best of our knowledge, our 
^^Cu NQR data in Fig. 1 is the first of its kind to de- 
tect the temperature dependence of the inhomogeneous 
distribution of the hole concentration in La2-a:Sra;Cu04 
or any other high Tc materials with quenched disorder. 

One important observation about Fig. 2a is that the 
B-linc ||] shows a very similar distribution in to 
the A-linc. The B-line itself originates from structural ef- 
fects |l^ . More specifically, a B-site corresponds to a Cu 
nucleus with a Sr^+ ion situated either directly above or 
below it, as supported in Fig. 2c where the fractional in- 
tensity of the B-line is shown to increase as x. We recall 
that the ^'^Cu NQR spectrum in La2-2;-j,Euj,Sra:Cu04 
for ^■^Cu isotope enriched samples |ll[ shows a third 
structural peak (the C-line) whose fractional intensity 
is equal to y. The C-line corresponds to a Cu nucleus 
directly above or below a Eu''+ ion. 



The similarity of the B-line to the A-line distribution of 
^^l/Ti rules out the possibilty that the distribution in the 
hole concentration Axhoie is concentric about a B-site. 
Rather, it indicates that the spatial modulation of Axhoie 
consists of patches that cover both A and B-sites equally 
and the size of each patch is larger than the average Sr^+ 
to Sr2+ distance of -1(2) nm for x = 0.16(0.04). 

We emphasize that the distribution in is not 

caused by flaws in our single phased ceramic samples. 
®^Cu NMR and NQR in high quality single crystals of 
both La2-2;Srj;Cu04 ||2^ and La2Cu04+5 |^ exhibit 
similar results. EMPA (electron micro-probe analysis) on 
single crystal Lai.85Sro.i5Cu04 showed that the spatial 
variation of the Sr^"*" content averaged over the focus area 
of the electron beam of ^ 1/im is Axsr — 0.01 <^ Axhoie- 
This rules out the possibility that /^m scale inhomo- 
geneities in the Sr^+ concentration is the cause of Ax hole- 
Moreover, earlier NQR measurements by Fujiyama et. al. 
|]lO| in La2-a;Sra;Cu04 without isotope enrichment agree 
with ours. We also note that the planar ^^O NMR Knight 
shift measurements provide additional evidence for a spa- 
tial distribution in the hole concentration. Our system- 
atic measurments for single crystals with x = 0.025, 
0.035, 0.05, 0.115 and 0.15 showed an overlap of the pla- 
nar -'^^O NMR central transition very similar to Fig. 2b. 
This means that the local static spin susceptibility has a 
similar spatial variation. 

We now turn our attention to the analysis of the 
^^Cu NQR spectrum. Experimental resolution limited 
the Axhoie data deduced from the distribution in 
above — 500K since all l/T}* curves start to merge above 
this temperature (Fig. 2). Not only does the spectral 
analysis enable us to extend Axhoie to higher tempera- 
tures, but we also get crucial insight into the length scale 
of the spatial variation Axhoie- 

The resonance frequency in NQR is proportional to the 
EEC surrounding the Cu nucleus (H] . The broad NQR 
linewidth shown in Fig. 2a and 2b implies that the lo- 
cal charge environment has a broad distribution in the 
EEC Possible causes of the broad distribution include: 
the random substitution of La'^^ with Sr^+, the spatial 
variation Axhoie , and the lowering of the local symmetry 
by lattice distortions. In order to simulate these effects, 
we carried out a point charge lattice summation to deter- 
mine the EEC at the Cu nuclear site. The conventional 
way of incorporating the Sr^+ substitution into an EEC 
point charge calculation is to reduce the effective charge 
of the La'^+ ion to La*^'^~^/^^+ However, in order to 
estimate the effect of the random replacement of La'^^ 
with Sr^"*", we ran a simulation of a realistic lattice with 
randomly distributed Sr^+ ions in a La'^+ matrix pi] . 
We found that the linewidth thus deduced is too narrow 
to account for the experimental result, as shown by the 
dashed line in Fig. 2a. 

Next, we ran the simulation with Axhoie 7^ using the 
following procedure: we first define an in-plane circle of 
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radius Rhoie surrounding the Cu nucleus a. The circle 
defines a "patch" of diameter 2Rhoie within which the 
local hole concentration is determined by counting the 
number TV of Sr^+ ions within the patch. The N donated 
holes from the TV Sr^+ ions are then uniformly distributed 
among the planar oxygen sites within the patch. The 
EFG contributions from all the point charges surrounding 
the Cu site a is then summed up, the corresponding NQR 
frequency deduced, and a count is placed in the appro- 
priate frequency bin. The lattice is then re-randomizcd, 
and the whole proceedure is repeated ^ 10'* times until a 
sufficient number of counts is present to reproduce the en- 
tire NQR spectrum, including the B-line. The aforemen- 
tioned case of the simulation with /\xhoie=^ (i-e. homo- 
geneous distribution of holes) corresponds to Rhoie ^ oo 
and yields a linewidth much narrower than the experi- 
mental data. Using Rhoie ~1.5 nm yields the dashed-dot 
spectrum in Fig. 2a, which is clearly too broad. Us- 
ing Rhoie ~3 nm however, results in the best fit shown 
by the solid lines in Fig. 2a and 2b. Once the best fit 
Rhoie is determined, Axhoie is directly obtained from the 
computation and is presented as the dashed lines in Fig. 
1. According to the simulation, the ^ 25 % increase in 
linewidth from 600K (Fig. 2b) to BOOK (Fig. 2a) for 
X = 0.115 can be accounted for by a ^ 25 % decrease in 
Rhoie which corresponds to 45 % increase in Axhoie 
from 600K down to 300K. We discuss implications of the 
temperature dependence later. 

We caution that we have ignored potential local tilt- 
ing of the CuOg octahedra caused by subtitution of the 
smaller Sr^+ ions, or potential precursor effects 25 - 27|] 
to the long range LTO (low temperature orthorhombic) 
transition. Since local octahedron tilting increases the 
NQR resonance frequency |2^Jl^, potential local distor- 
tions of the lattice may provide the additional mecha- 
nism to broaden the spectrum. We have ignored this 
and attributed the extra line broadening to Axhoie, 
hence Axhoie thus deduced may be an upper bound 
and Rhoie ~ 3 nm maybe a lower bound. We note 
that the minimum patch size of 2Rhoie >6 nm is still 
larger than the average Sr^+ ion separation of 1(2) nm 
for X = 0.16(0.04), which further rules out the possibilty 
that the distribution of the local hole concentration is 
concentric about each Sr^+ ion. The validity of our sim- 
ulation breaks down altogether in the temperature region 
where local octahedron tilting actually becomes apparent 
due to long range LTO transition, as indicated by an in- 
crease in ^^VQ ||l^,|l^. The analysis shown in Fig. 1 is 
thus terminated below this temperature. 

Our lineshape analysis sets only the lower bound of 
the length scale 2Rhoie >6 nm of the spatial variation 
Axhoie — 0.03 ^ 0.06. On the other hand, our afor- 
mentioned EMPA analysis indicates that the deviation 
in Sr^+ concentration is as small as Axsr ^ 0.01 aver- 
aged over the length scale of ^ 1/im. Furthermore, ear- 
lier high resolution X-ray diffraction studies p3 showed 



that Axsr ^ 0.01 averaged over the length scale of > 10 
nm. These results imply that the local hole concentra- 
tion Xhoie modulates as much as Axhoie = 0.03 ~ 0.06 
within the short distance scales of 6-10 nm, but the spa- 
tial average of Xhoie over greater length scales is as little 
as ~ 0.01. One obvious possibility is the clustering of 
Sr^+ ions with very short length scale of 6-lOnm. How- 
ever the absence of any substantial diffuse scattering in 
neutron diffraction experiments | |2^ makes this scenario 
unlikely. Moreover, the presence of the temperature de- 
pendence of Axhoie suggests that the temperature inde- 
pendent quenched disorder caused by the clustering of 
the Sr^+ ions alone cannot entirely account for our re- 
sults. We recall that our lineshape analysis is based on 
a completely random distribution of Sr^+ ions. We also 
note that the spin-spin correlation length ^ cannot be 
identified as Rhoie because ^ is much shorter than 2Rhoie 
and increases with decreasing temperature ]30|] . The lat- 
ter qualitatively contradicts with the decrease of Rhoie 
with decreasing temperature. These considerations in- 
evitably lead us to conclude that there must exist an 
electronic mechanism at short length scale which causes 
the holes to segragate. 

One possible mechanism of establishing such a short 
electronic length scale is phase separation (which may or 
may not have a stripe- like structure). In fact, Ino et. al. 
pl[ reported evidence for phase separation at ~ 80K in 
La2_a;Sra;Cu04 for x <0.15 based on photoemission ex- 
periments, i.e. the chemical potential does not vary with 
increasing hole doping up to x = 0.15. In this concentra- 
tion regime, the presence of a stripe instability has been 
firmly established at low temperatures |^,^,|ll|. Our 
new observation in Fig.l may be a precursor to such a 
phenomena. However, we observed comparable or larger 
values of Axhoie for x — 0.20 and 0.25 where photoemis- 
sion measurements ||3l|] do not detect signatures of phase 
separation. Thus an additional mechanism(s) for segre- 
gation of holes is necessary to account for our results at 
least above x ~ 0.15. Interestingly, the temperature de- 
pendence of Axhoie and the magnitude of Rhoie may be 
qualitatively understood based on a toy model of a ther- 
mal activation type process for the segragation of holes. 
Using the in-plane static dielectric constant of ~ 30 
and the effective mass m* 2 (in units of TOe) deduced 
from lightly doped La2Cu04+5 132], one obtains a factor 
m* /el ^ 1/450 smaller binding energy of E;, ~ 30meV 
or Eh ^ 350K and a factor €s/m* ~ 15 larger Bohr ra- 
dius of ao ^ Inm. These temperature and length scales 
are indeed comparable to our observation. Even though 
this oversimplified model docs not include the effects of 
screening for instance, it might suggest that the donated 
holes are bound by a Coulomb type potential created by 
the Sr^^ ions rather than homogeneously doped into the 
Cu02 planes. Needless to say, our toy model analysis 
employed here above '--^lOO K does not necessarily sup- 
port or rule out the presence of intrinsic phase separation 
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at lower temperatures. 

The work was supported by NSF DMR 98-08941 and 
99-71264. 



[32] C.Y. Chen et. al., Phys. Rev. B 43, 392 (1991). 

FIG. 1. Temperature dependence of the distribution in lo- 
cal hole concentration /\Xhoie as deduced by ®^Cu NQR l/Tf*^ 
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FIG. 2. (a) Frequency dependence of "^^l/Ti (□) across 
•^^Cu NQR spectrum (•) in Lai.885Sro.ii5Cu04 at 300K. 
Dashed, solid and dashed-dot spectra are theoretical fits with 
patch sizes Rhoie ^ 00, 3nm and 1.5nm respectively. Dotted 
lines through "^^l/Ti are a guide for the eye. (b) '^^Cu NQR 
spectra in a; = 0.04 (A), x = 0.07 (o), x = 0.115 (•), and 
X — 0.16 (a) at 600K. Solid lines are theoretical fits with 
Rhoie ~ 3nm. (c) Fraction fs of B-line intensity to total in- 
tensity (with T2 corrections) at (+) 600K and (x) 300K as a 
function of x. 



FIG. 3. Temperature dependence of l/Tj (x), 1/T~^ 
(▼), 1/rf (o), and 1/T+ (A) in La2-^Sr^Cu04 where x is 
given in each panel. Solid lines are a guide for the eye, and 
dashed lines are for x = 0.00, 0,02, 0.04, 0.07 ,0.09, 

0.115, and 0.16 where 1/Tf monotonically decreases with in- 
creasing X. 
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